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ABSTRACT 

We establish that the Magellanic Stream (MS) is some 40° longer than previously known with 
certainty and that the entire MS and Leading Arm (LA) system is thus at least 200° long. With the 
Green Bank Telescope, we conducted a ~200 deg 2 , 21-cm survey at the tip of the MS to substantiate 
the continuity of the MS between the Hulsbosch & Wakker data and the MS-like emission reported 
by Braun & Thilker. Our survey, in combination with the Arecibo survey by Stanimirovic et al., 
shows that the MS gas is continuous in this region and that the MS is at least ~140° long. The 
MS-tip is composed of a multitude of forks and filaments. We identify a new filament on the eastern 
side of the MS that significantly deviates from the equator of the MS coordinate system for more 
than ~45°. Additionally, we find a previously unknown velocity inflection in the MS-tip near MS 
longitude Lms~ —120° at which the velocity reaches a minimum and then starts to increase. We 
find that five compact high velocity clouds cataloged by de Heij et al. as well as Wright's Cloud are 
plausibly associated with the MS because they match the MS in position and velocity. The mass of 
the newly-confirmed ~40° extension of the MS-tip is ~2 x 10 T M (d/120 kpc) 2 (including Wright's 
Cloud increases this by ~50%) and increases the total mass of the MS by ~4%. However, projected 
model distances of the MS at the tip are generally quite large and, if true, indicate that the mass of 
the extension might be as large as ^10 8 M©. From our combined map of the entire MS, we find that 
the total column density (integrated transverse to the MS) drops markedly along the MS and follows 
an exponential decline with Lms of 7Vhi=5.9 x 10 21 exp(LMs/19.3°) cm -2 . Under the assumption 
that the observed sinusoidal velocity pattern of the LMC filament of the MS is due to the origin of 
the MS from a rotating LMC, we estimate that the age of the ~140°-long MS is ~2.5 Gyr. This 
coincides with bursts of star formation in the Magellanic Clouds and a possible close encounter of 
these two galaxies with each other that could have triggered the formation of the MS. These newly 
observed characteristics of the MS offer additional constraints for MS simulations. In the Appendix 
we describe a previously little discussed problem with a standing wave pattern in GBT Hi data and 
detail a method for removing it. 

Subject headings: Galaxies: interactions - Galaxies: kinematics and dynamics - Galaxies: Local Group 
- Galaxy: halo - Intergalactic Medium - Magellanic Clouds - ISM: Hi 



1. INTRODUCTION 

More than three decades ago, evidence of a galactic 
Hi stream around the Milky Way (MW) was discovered 
as an ~100°-long Hi tail emanating from the Magel- 
lanic Clouds (MCs) and named the Magellanic Stream 
(MS; lWannier fc Wrixon||1972| |Mathewson et al.j|1974| ). 
The Leading Arm (LA), a counterpart to the MS that 
stretches in front of the MCs, has also been known for 



many decades and suspected to be from the MCs ([Math 

;l y linki 

the MCs until t he studies by Put man et alj ([1998 ) and 



ewson et al. 1974), but was not definitively linked to 



Lu et al. 



structure 



(1998). The complex and double-fila mentary 
of the MS was fir st clearly shown by Wakker 



([2001 ), based on the data of |Hulsbosch fc Wakker[ ( 1988 
hereafter HW88), and later sho wn in more detail in the 



higher-resolution HIP ASS data ( |Barnes et al.||2001| ) by 
|Putman et all ([20031 hereafter P03 )~ 

Recently, jBraun & Thilker| ( |2004[ hereafter BT04) con- 
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ducted a deep and wide-field Hi survey of the M31 region 
with the Westerbork Synthesis Radio Telescope (WSRT). 
They serendipitously discovered some faint emission in 
the western portion of their survey that appeared to be 
consistent with a spatial and kinematical extension of the 
MS, but at much lower column densities (Fig. [I]). This 
suggested that the MS was plausibly longer than previ- 
ously recognized, although a ~10° gap existed between 
the MS as revealed by P03 and the faint emission discov- 
ered by BT04 (Fig. [l] as well as Fig. 5 of BT04). While 
the HW88 data show the MS continuing to the edge of 
the BT04 survey, these data lacked sufficient resolution 
and coverage to be certain (Fig. [I]). A better map of this 
region would help to clarify the continuity of the MS in 
this region. 

Recent simulations by |Besla et al.| ( |2007|) , mo tivated by 
new prope r motion measurements for the MCs ( |Kallivay^] 
alii et al.|2QQ6a|b| |Piatek et al.|2008[ ), conclude that the 
MCs might be on their first passage around the MW. 
The new orbit predictions suggest that the MCs have 
spent most of their lives isolated from the MW, a sce- 
nario that has major consequences for the viability of the 
two primary MS mechanisms that have traditionally been 
postulated for the formation of the MS - ram pressure 



stripping (e.g., |Meurer et al.||1985l |Sofue|1994 
Davis|1994}|Mastropietro et al.|2005| ) and tidal 



Moore 



stripping 
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Fig. 1. — HI detections of the Magellanic Stream in the region 
of our survey by previous authors: open blue circles indicate MS 
detections by HW88; the grayscale image shows TVhi from P03 with 
contour levels of 0.2, 0.5, lxlO 19 atoms cm 2 ; and the pink image 
shows TVhi from BT04 with contour levels at 0.2, 1, 5xl0 18 atoms 
cm 2 . The dotted line shows the edge of the P03 data; the dashed 
line, the edge of the BT04 data. 



Yoshizawa & Noguchi 2003; Connors et al.|2006[). To re- 



Murai fc Fujimoto|1980l|Gardiner fc Noguc hi 1996 



move the MS gas effectively, both of these mechanisms 
require the MCs to be relatively close to th e MW for a 
prolon ged encounter time. If, as suggested by |Besla et al.| 
( |2007| ), the orbits of the MCs have kept them more iso- 
lated in the past, then it becomes difficult for either ram 
pressure or tidal forces alone to create the MS. The find- 
ing of an even longer MS would exacerbate this problem 
for the two standard mechanisms since a greater length 
would require gas to be pulled out of the MCs at even 
larger MW-distances where the density of a ram pressure 
medium an d tidal forces are eve n lower. 

However, |Nidever et al. ( 2008[ ) put forward a third MS 
formation mechanism after (1) tracing one filament of the 
MS as well as the LA back to their origin in the southeast 
Hi overdensity (SEHO) in the LMC, and (2) illustrat- 
ing that the SEHO contains large gaseous outflows from 
supergiant shells. The "blowout hypothesis" postulates 
that supergiant shells in the dense SEHO blow out gas 
from the LMC to large radii where it is easier for ram 
pressure and/or tidal forces to fully strip the gas and 
disperse it to create the MS and LA. Because blowout 
does the hard work of overcoming a significant portion 
of the gravitational hold of the LMC on the gas, this sce- 
nario should be an effec tive solution to the l arge-distance 



dilemma raised by the [Besla et al. (2007) fir st-passage 
scena rio. I n a proposed alternative scenario, |Mastropr 
( |2009| ) was able to produce a 120°-long MS with 
ram pressure (and tidal forces) using a high- velocity, un- 
bound LMC orbit, but whether this new model can ac- 
count for an even longer MS remains unclear. A better 
understanding of the true length of the MS would place 
important, more specific constraints on these various MS 
formation scenarios. 

In this paper we address the question of whether the 
BT04 emission really is an extension of the MS and if 
the MS is actually much longer than previously widely 



accepted. In pursuit of this goal we conducted a 200 
deg 2 Green Bank Telescope (GBT; |Lockman|[T998 ) 21- 
cm survey in the region between the tip of the classical" 
MS and the features found in the BT04 survey to explore 
whether there is continuity of the MS across this region. 
O ur GBT data, in combinat ion with the Arecibo survey 
of Stanimirovic et aT| ( |2008| hereafter S08), demonstrate 
that the MS is both spatially and kinematically contin- 
uous across this region and therefore that the MS is at 
least 140° long (~40° longer than previously verified), 
and probably even longer. This additional extension of 
the MS makes the entire MS system (including the LA) 
at least 200° long. In addition, we combine all available 
Hi data of the MS-tip into one datacube to investigate 
the structural characteristics of the MS in this region. 
The MS-tip is composed of many thin filaments, as pre- 
viously noted by S08, but we identify a new filament in 
the eastern portion of the MS that markedly deviates 
from the equator of the MS coordinate system for some 
45°. Additionally, we find that there is a previously un- 
known velocity inflection at the MS-tip where the MS 
radial velocities reach a minimum and then begins to in- 
crease. The mass of the newly mapped portion of the 
MS-tip is -2 x 10 7 M (d/120 kpc) 2 which increases the 
estimated mass of the entire MS by ~4%. 

This paper is organized as follows: A brief descrip- 
tion of our GBT survey and data reduction is given in 
[J2j while more details on the data reduction and some 
difficulties we encountered in dealing with GBT 21-cm 
spectral data are presented in the Appendix. In ^3] we 
combine the various extant Hi datacubes of the MS tip 
and analyze the structural characteristics of the extended 
MS. Figures and descriptions of the entire 200° MS+LA 
system are presented in Q A discussion of our findings 
is given in §5\ and a summary of our primary conclusions 
is given in q6| 

2. GBT SURVEY 
2.1. Observations 

We used the GBT to conduct an ~200 deg 2 , 21-cm sur- 
vey (proposal ID: GBT06A-066, 102 hours) to study the 
continuity of the MS in the region between the "classical" 
MS and the BT04 survey. The survey area is centered 
on (l,b) « (94°, -40°) and extends (~28°x22°) in (I x b) 
(but not completely filled; see Fig. |2j. Our survey area 
was partly chosen to be complementary to the Arecibo 
survey of Stanimirovic et al. We used the "On-the-Fly" 
mapping mode (scanning in Galactic longitude) to ob- 
tain frequency- switched, 21-cm spectral line data with 
the Spectrometer back end. The usable velocity range 
is —960 < Vlsr < +540 km s _1 , with a velocity resolu- 
tion of 0.16 km s _1 . The observations were taken in 49 
"bricks" of dimension 2.1° x 2.1°. Each brick consists of 
a 37x37 array of 4-second integrations, with a 3.5 ; spac- 
ing between each integration (which fully samples the 
GBT 21-cm half-power beam width of 9.2 / ). Neighbor- 
ing bricks were overlapped by 6' to allow for a consistent 
calibration across the survey area. The survey area can 
be seen in Figure [2] 

During each observing session calibration data were 
obtained for th e IAU standard positions S7 and S8 
(Williams 1973[) , and, for our own secondary standard 
position, at = (103.0°, -40.0°). 
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Fig. 2. — (a) HI column density, 7V H i, (in units of 10 18 atoms cm -2 ) of the MS gas (-502 < V LS r < -266 km s -1 ) after 3cr filtering of 
the GBT datacube. The end of the detected MS-tip from the P03 HIPASS data is shown in white contours (at 0.2, 0.5 and lxlO 19 atoms 
cm -2 ), while gray contours show the MS-like gas in the southern edge of the BT04 survey (at 0.2, 1, and 5xl0 18 atoms cm -2 ). There 
is nearly-continuous MS emission across our survey at I « 93° and I « 100°. This shows continuity of t he M S in the eastern part of the 
region (between the "classical" MS and the BT04 survey) where continuity was previously uncertain (see £ |3.2| ). The northwestern region of 
our survey is filled with MS cloudlets and suggests that the western part of the MS is also continuous across this region, (b) First-moment 
map for the data in (a). Hue indicates the flux- weighted velocity (Vlsr m km s — 1 as indicated by the colorbar); intensity indicates Nm. 
The velocity of the I « 100° gas is ~70 km s -1 lower than the gas at I w 93°. A dashed line demarcates the boundary of our survey. The 
figures have been oriented so that they are aligned with the other MS figures in this paper (e.g., Figs. ^] and |5|. 



2.2. Data Reduction 

We used the GETFS program in GBTIDL 5 to calibrate 
our frequency- switched data using a basic T a = T^l x 
(sig — ref)/ref calibration. Each integration and linear 
polarization was calibrated separately. 

We used our own special-purpose IDL routines to per- 
form the baseline fitting and removal. Before baseline 
removal, the data were binned 10 x (giving a velocity 
resolution of ~1.6 km s _1 ) to decrease the noise. The 
baseline was fitted for each integration and polarization 
separately. For the XX polarization, a 5th-order polyno- 
mial was fitted to the 21-cm spectrum but excluding the 
Galactic emission at —140 < Vlsr < +100 km s -1 . The 
fitting process was done iteratively so that emission lines 
(anything above 2<j in a Gaussian smoothed version of 
the spectrum) could be excluded from the fit. The pro- 
cedure for the YY polarization was similar although it 
included a sinusoidal component in order to remove a 
standing wave with a period of ~1.6 M Hz in the GBT 
data (this is explained in more detail in j |A.l| ). 

After baseline removal, the median (calculated from 
an emission-free region of the spectrum) was subtracted 
from each polarization and then the two polarizations 
were averaged to obtain the final, reduced spectrum for 

5 http://gbtidl.nrao.edu/ 



that position. The RMS noise is ~0.072 K per 1.6 
km s _1 channel; our 3<r sensitivity is ~4.9xl0 18 cm -2 
over 20 km s -1 . 

All 49 bricks were regridded onto a single datacube 
with a Galactic cartesian grid and a 3.5' spacing using the 
IDL routine TRIGRID with linear interpolation (each ve- 
locity channel independently). The datacube was then 
Gaussian smoothed in velocity (FWHM = 16 km s _1 ) 
and spatially (FWHM = 10.5'), giving a final RMS noise 
level of ~0.0075 K per 1.6 km s _1 channel. Our 3a sen- 
sitivity is ~5.2xl0 1T cm -2 over 20 km s _1 in the final, 
smoothed datacube. 

We found and removed a persistent but weak (~5mK) 
emission line in our datacube. We call this feature the 
GBT Hi "ghost" and find that it is spurious and likely 
due to an improperly filtered sideband on a local oscilla- 
torat the GBT. A more detailed description is given in 



£A.2 



2.3. GBT Survey Results 

Figures [2] and [3] show our final, smoothed GBT dat- 
acube. An "RMS map" that contains a robust measure 
of the RMS noise for each position was also created (not 
shown) . To pull out real features the data were a-filtered 
so that only pixels in the datacube that were above 3x 
the RMS noise (for their respective position) were used 
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Fig. 3. — (a) Integrated intensity of the datacube along b (in units of K deg) after cr-filtering. The velocity of the I « 100° gas 
(Vlsr ~ —406 km s _1 ) is ~70 km s _1 lower than the gas at I w 93° (Vlsr ~ —337 km s _1 ). (b) Integrated intensity of the datacube 
along I (in units of K deg) after cr-filtering. The velocity gradient of the MS with b is clearly apparent (since b is nearly parallel with the 
MS at this location). Some of the I w 100° cloudlets can be seen slightly below the main linear trend of MS gas at b ~ —49° and —40°. 



to make the figures. 

The column density of the MS gas (-502 < Vlsr < 
—266 km s _1 ) is shown in Figure j2a. Contours of the 
end of the MS from the P03 HIPASS data (white) and 
the MS-like gas in the southern edge of the BT04 survey 
(gray) reveal the region where the continuity of the MS 
is unclear and that is filled by our GBT survey (mostly 
on the eastern part) 6 . The HW88 data detect the MS in 
a wedge shaped region in the central part of this figure 
(Fig. IT]). There is nearly-continuous MS emission across 
our G13T survey at I ~ 93° which confirms the conti- 
nuity of the MS in this region, as first seen by HW88, 
and reveals it to be composed of many small clumps and 
cloudlets. Additionally, there is nearly-continuous MS 
emission across our survey at I « 100° that connects 
the "classical" MS to the BT04 survey for the eastern 
(I « 100°) portion of the MS. Previous data (HW88, 
P03) did not show continuity of the MS in this east- 
ern region. The northwestern region of our survey is 
filled with cloudlets at MS velocities and suggests that 
the western portion is also continuous across this region. 
Figure shows the velocity map for the MS gas. The 
gas at I « 100° is at slightly lower velocity (more neg- 
ative) than the gas at I « 93°. Figure shows the 
integrated intensity of the datacube along and Figure 
[3)3 the integrated intensity along I. The velocity gradient 
of the MS with b is clearly apparent (since the MS nearly 
follows a constant line of Galactic longitude at this loca- 
tion). Some of the I ~ 100° cloudlets can be seen slightly 
below the main linear trend of MS gas at b ~ —49° and 
-40°. 

In summary, our dataset shows that the MS is continu- 
ous across this region between the end of the "classical" 
MS and the BT04 survey. Therefore, the MS is much 
longer than previously known with confidence. In the 

6 In this region of the sky the celestial and galactic coordinate 
systems are nearly aligned and, therefore, references such as "east- 
ern" or "northwestern" can refer to either coordinate system. 



next section we combine our data with other Hi datasets 
in order to obtain a complete picture of the tip of the 
MS and to determine the full extent and distribution of 
the MS on the sky and in velocity. 



3. 



THE COMBINED MAGELLANIC STREAM TIP 
DATACUBE 



To gain a clear picture of the entire tip of the MS we 
combined our GBT datacube with (1) the S08 Arecibo 

)5, hereafter 



datacube, (2) a part of the Briins et al. d2 
Br05) Parkes datacube, and (3) the BTU4 Westerbork 
datacube. 



3.1. Combining the GBT, Arecibo, Parkes, and 
Westerbork Datacubes 

Our combined MS-tip datacub e is on a Cartes i an gri d 
on the MS coordinate system of Nidever et al.| ( |2QQ8| ) 7 
with an angular step size of 3', a velocity step of 1.5 
km s _1 , and a velocity range of —548.5 < Vlsr < —100.0 
km s _1 . A cubic spline was used to interpolate the dat- 
acubes onto the new velocity scale and the IDL function 
TRIGRID was used to regrid the datacube spatially one 
velocity channel at a time. Each dataset was regridded 
onto the new MS-tip grid separately before all datasets 
were combined. 

GBT Data: The raw GBT data have an angular resolu- 
tion of 9.2', angular sampling of 3.5', and velocity spacing 
of 0.16 km s _1 (before binning). The original, velocity- 
binned, GBT data were griclded onto the new dat- 
acube, and then smoothed with a Gaussian (FWHM=15 
km s _1 ) in velocity and spatially (FWHM=6 / ). The av- 
erage final RMS noise is 0.011 K. 

Arecibo Data: The S08 Arecibo datacube is a combina- 
tion of three observing programs with differing S/N. The 

7 Note that this coordinate system was designed to bisect the 
Magellanic St ream and differ s from the "Magellanic coordinate sys- 
tem" used by Wakker] ( pOOl) and P03. 
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original data have an angular resolution of 3.5', angular 
sampling of 1.8', and velocity spacing of 0.18 km s _1 , but 
after gridding have an angular spacing of 3' and velocity 
spacing of 1.47 km s _1 (after binning). We masked out 
regions of the datacube with missing or bad data, as well 
as regions with poor S/N that overlapped the GBT data 
(mainly in the northern part). After the Arecibo dat- 
acube was regridded onto our MS-tip datacube it was 
smoothed with a Gaussian (FWHM=15 km s _1 ) in ve- 
locity. The average final RMS noise is 0.012 K. 

Parkes Data: The Br05 Parkes data have an angular 
resolution of 14.1, angular sampling of —5', velocity 
spacing of 0.82 km s~\ and cover a velocity range of 
—410 < Vlsr ^ 380 km s _1 , but after gridding have an 
angular step of 10' with the same velocity spacing. Af- 
ter the Parkes datacube was regridded onto our MS-tip 
datacube it was smoothed with a Gaussian (FWHM=15 
km s _1 ) in velocity. The average final RMS noise is 0.032 
K. 

Westerbork Data: The BT04 Westerbork auto- 
correlation data have an angular resolution of —50' (ef- 
fective telescope beam), an angular sampling of 15', and 
velocity spacing of 8.3 km s _1 (Hanning smoothed to 
give a spectral resolution of 16.6 km s _1 ). The gridded 
datacube was resampled onto the Local Group Standard 
of Rest (LGSR) velocity system with a velocity range 
of -535 < Vlqsr < 95 km s" 1 (-700 < V LS r < -50 
km s _1 ). After converting the velocities to Vlsr, we dis- 
covered that there was a systematic velocity offset be- 
tween the Westerbork data and the GBT and Arecibo 
data where they overlapped. We use d the Leiden- 



Argentine-Bonn (LAB) Hi all-sky survey ( jKalberla et al. 
2005 ) to determine the velocity offset in all positions that 
had emission detected in both the LAB and Westerbork 
data. The velocity offset changed systematically with po- 
sition with an average offset of —30 km s _1 . The additive 

Offset (i.e. Correct = ^Westerbork + ^offset) Was well-fitted 

by 

Voffset = 58.24 + 0.1954a - 0.00959a 2 - 0.74355 

where a and S are in degrees. We applied this velocity off- 
set to the Westerbork data and converted to Vlsr before 
regridding onto our MS-tip datacube. Since M31 and 
M33 are very prominent in the Westerbork datacube and 
would "contaminate" our MS position-velocity diagrams, 
we blanked their emission from the datacube. The aver- 
age final RMS noise is 0.0012 K. 

We used a rank order of decreasing priority as 
GBT:Arecibo:Parkes:Westerbork to select which data to 
use for regions where the datasets overlap. We cre- 
ated a "hi-res" datacube, where all four datasets have 
their original resolutions (Fig. Hk), and a "low-res" dat- 
acube, where the combined GBT /Arecibo/ Parkes data 
were smoothed to a resolution of 50' before being com- 
bined with the Westerbork data (Fig. [4)3). In the hi-res 
datacube the intensity of the emission in the Wester- 
bork part is generally low (by — 2x) due to beam dilu- 
tion, which creates an artificial drop in emission/column 
density at the interface between the Westerbork and 
GBT/ Arecibo data. The low-res datacube gives gen- 
erally smoother images because the data are relatively 
homogeneous in angular resolution. 

We have also identified five Compact or extended High 



Vel ocity C louds (C/HVCs) from the catalog of |de Heij et 
aDd2062) (numbers 237, 304, 307, 391, 402 in their Table 
TJtnat lie outside the area covered by the four datasets in 
our MS-tip datacube but might be associated with the 
MS because they match the MS in position and veloc- 
ity. The C/HVCs were added into the combined MS-tip 
datacubes with the structural Hi parameters given by de 
Heij et al. Since the Leide n/Dwingeloo Hi Survey (LDS; 
Hartmann fc Burton|1997|) was used to identify these Hi 



clouds, they have an angular resolution of 36' and ve- 
locity resolution of ~1 km s _1 . Most of the CHVCs of 



Westmeier fc Koribalski] fl2008| hereafter WK08), which 
they identified as being associated with the MS, also lie 
outside the area surveyed by the four datasets. However, 
structural parameters are not available for these CHVCs 
and therefore we did not add them to our MS-tip dat- 
acube but instead plot them as dots in Figures [4}T~ 



3.2. MS Tip Results 

The Hi column density map of the hi-res/low-res data 
can be seen in Figure [4] after an MS velocity selection 
(removing the well-separated local Galactic and interme- 
diate velocity gas - see Fig. |6k) and 3<r filter have been 
applied. With the addition oithe GBT and Arecibo data 
it is now clear that the entire MS is continuous with the 
emission identified by BT04 as an extension of the MS 
(starting at 5=+20°). The MS is at least -140° long 
and thus is ~40° longer than previously established with 
certainty. The entire MS system — including the LA 
— is therefore at least ~200° long. The MS extends to 
the limit of the MS-tip datacube (the Westerbork data) 
at S ~ +50° and therefore is plausibly even longer than 
shown here. In fact, there is even one HVC from |de Heij| 
let all (120021) (number 424) that is consistent in position 
and velocity with being an extension of a filament on the 
eastern side of the MS at Lms~ —163° (see Fig.[7]below). 

The Hi data of the MS-tip (Fig. H) show a complex 
richness of patterns that look like forKs and multiple fil- 
aments. The MS splits and diverges into two filaments 
near Lms ~ —80°, which were already visible in previ- 
ous data (HW88, P03), but there are more forks farther 
along the Stream. S08 noted that the western region of 
the MS-tip separates into four thin filaments that they 
named S1-S4. The four S08 filaments are indicated in 
Figure [4^i and are extended slightly now that the BT04 
survey data are included. S3 and S4 appear to run off the 
edge of the survey coverage, while SI and S2 are likely 
converging or overlapping for Lms ^ —108°. 

A new filament of the MS (on the eastern side) is visi- 
ble in our data that separates from SI near Lms ~ —95°. 
Following the S08 nomenclature we call this filament SO. 
It is made up of many cloudlets, especially right after 
it separates from SI, but the continuity of SO is clear 
in our GBT data (Fig. [2|. As can be seen in Figure 
[4j SO extends all of the way to the end of the cover- 
age area near Lms~ —140°. It deviates from the equa- 
tor of the MS coordinate system and extends diagonally 
at an angle of -14.6° (slope of -0.26). The SO fila- 
ment also follows a great circle with a north pole of 
(Z,6)=(205.2° -15.4°) (Fig. |7j black-white dashed line) 
which is —18° from the nortn pole of the MS coordinate 
system (/,6)=(188.5°,— 7.5°). The western portion of the 
MS, which is probably a combination of SI, S2, and pos- 
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Fig. 4. — HI column density, A^i, of the MS-tip combined datacube (log(A^Hi) in units of atoms cm -2 ) after an MS velocity selection 
and 3cr filtering have been applied, (a) The "hi-res" datacube with all four datasets (GBT, A recib o, Parkes, and Westerbork) in their 
native resolutions. A celestial coordinate grid is overlaid. The various MS filaments described in §3.2| are outlined with dashed lines. S1-S4 
were defined in S08 and we extend them slightly with the addition of the BT04 data. We define the new eastern-most filament as SO; 
which splits with SI near Lms~ —95° but extends until Lms~ —143°. (b) The "low-res" datacube with all datasets at 50' resolution. 
A Galactic coordinate grid is overlaid and the boundary of our GBT survey is indicated by a solid white line. CHVCs from WK08 that 
they associated with the MS are represented by white dots. A dashed line indicates the boundary of the combined HI datacubes, and a 
dashed-dotted line delineates the eastern edge of the region studied by WK08 (0 h 00 m < a < l h 30 m and -30° < S < +25°). The MS 
is seen here to be ~40° longer than the "classical" MS (which ends near Lms ~ — 105°), and now extends to the end of the Westerbork 
survey at 5 « +50°. The MS is very complex in this region with many forks and filaments — see text for details. Wright's Cloud is the 
backward "L"-shaped cloud at (L M s^ms)~(-114 ,+29°). 
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Fig. 5. — Velocity map of the SO'-smoothed MS tip datacube. 
Hue indicates the flux weighted velocity (Vlsr m km s _1 as in- 
dicated in the colorbar) and intensity indicates Nm. The colored 
dots are CHVCs from WK08. The velocity gradient along the MS 
is clear from the color coding. 



sibly S3 8 , also extends to the end of the coverage area 
near Lms~ —140°. More data are needed in the west 
and north to track the MS filaments further. 

Figure [5] shows the velocity map of the low-res dat- 
acube with the WK08 CHVCs plotted as colored dots. 
The strong velocity gradient (~6 km s _1 deg~^) with 
Lms is apparent from the color coding. Figure [6] shows 
various Vlsr^ms position- velocity diagrams of the MS- 
tip low-res datacube. Figure[6^ shows the total intensity 
of the entire datacube integrated in Bms- The strong 
velocity gradient with Lms is again apparent, and so is a 
velocity inflection near Lms~ —120° where the velocity 
of the MS levels off and then starts to increase. This 
is the first evidence of a leveling-off of the MS velocity 
anywhere along its length. 

Figure [SJd is a Vlsr-^ms diagram similar to Figure 
[6^i, but only for the western portion of the datacube 
T^ms ^ 0°) showing the velocity behavior of the S2- 
S4 filaments and part of SI. Most of the gas follows the 
velocity inflection (and the velocity fiducial line) includ- 
ing the low column density gas at — 132° < Lms ^ —117°. 
There is a cloud at (L M s,^lsr)-(-135°,-480 km s" 1 ) 
that is consistent with a near-linear extrapolation of the 
MS, and that might indicate that not the entire MS 
participates in the velocity inflection. The clouds at 
(£ms,Vlsr)~(-120°,-320 km s" 1 ), above the main MS 
in the figure, are in the far western part of the datacube 
at Bms ~ —20° and are not part of the main MS filaments 
(S0-S4). 

A Vlsr~£ms diagram of the central portion of the dat- 
acube (0° < B MS < +15°) including filament SO (and 
parts of SI) is shown in Figure [6^. The velocity levels 
out and then increases, but not as quickly as for the west- 
ern portion of the MS (Fig. [6b). The velocity scatter for 
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S4 runs off the coverage near (£ms^ms)=( — 116°, — 11°). 



Fig. 6. — (a) Total intensity of the SO'-smoothed low-res MS-tip 
datacube integrated along Bms ( m units of K deg on a logarithmic 
scale). There is a velocity inflection near Lms ~ —120° where 
the velocity flattens out and begins to increase. White dots are 
the CHVCs from WK08, and a dashed fiducial line (same in all 
panels) shows the proposed general velocity trend. Wright's Cloud 
is the red blob at (L M s,Vlsr)~(-114°,-380 km s" 1 ). (b) Same 
as (a) but only for the western portion of the MS-tip datacube 
(Bms S 0°) showing the S2-S4 filaments and part of SI. (c) Same 
as (a) but only for the central portion of the MS-tip datacube (0° < 
Bms 5 +15°) sh owing filament SO (and parts of SI). Two of the|de| 
|Heij et al | (2002) CHVCs are at L M s ~ -140°. (of) Same as (a) but 
only tor the eastern portion of the MS-tip datacube (Bms > +15°) 
showing the clouds around Wright's Cloud and many of the WK08 
CHVCs. The sudden cutoff in velocity at Vlsr= — 410 km s _1 for 
Lms 92° in all panels is due to the velocity limit of the Br05 
Parkes datacube. 

the central region is larger than for the western. There 
are two clouds near (I/ms 5 Vgsr)~( — 140°,— 450 km s _1 ) 
that might not participate in the velocity inflection but 
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MS indicate that the distance increases very rapi dly at 




the tip of the MS especially for Lms ^ , — 100° (e.g., Gar- 



0.001 



0.01 



Fig. 7. — (a) Same as Fig. |4p but with an exten ded 
Lms range and showing HVC 424 from |de Heij et al.l ( |2002| ) at 
Lms~~ 163° that might belong to an extension ot the bO nTament 
of the MS. The black-white dashed line indicates the great circle 
with north pole of (Z,&)=(205.20°,-15.40°) that best fits the SO 
filament of the MS and which also passes close to HVC 424. The 
great circle also fits the eastern portion of the MS (S0+S1 before 
they split) all the way to Lms ~ —80°. (b) A similarly extended 
version of Fig. l6t. A second Vlsr^ms fiducial is shown that fits 
HVC 424 . 

continue to more negative velocities. Figure [7] extends 
the Lms range of Figure [B^ to include HVC 424 from [de| 
Heij et al. (2002) which appears to be an extension of 
the SO filament at Lms~ —163°. This cloud, as well as 
the clouds at Lms~ —140°, are consistent with a possible 
shallower velocity inflection for the SO filament than for 
the rest of the MS. 

Finally, Figure |6]i shows the part of the datacube to 
the east of the S O filament (L> ms ^ +15°) which includes 
Wright's Cloud ([Wright |1979D at L M s ~ -114° and most 
of the WK08 CHVCs. The clouds around Wright's Cloud 
appear to be an extension of the WK08 CHVCs and they 
all generally follow the MS velocity trend. Some authors 
(Wright 1979, BT04, Putman et al. 2009) have suggested 
that Wright's Cloud might be associated with the MS; 
the velocity agreement of the MS and Wright's Cloud 
in Figures [6^i and d makes this association likely. We 
discuss this possibility further in ^5] 

To calculate the mass of the newly-confirmed part of 
the MS-tip we remove all of the Hi emission that was pre- 
viously observed (i.e., the "classical" MS) by subtracting 
the Br05 Parkes data from our combined MS-tip column 
density map. Orbits of the MCs and simulations of the 



diner & Noguchi 199 6jjConnors et aL||200gj |Besla et al 



2007; Mastropietro|2009|), and the newer orbits and sim- 



ulat io ns" based on the new HST proper motions of the 
MCs ( jKallivayalil et a l. 2006a b) give significantly larger 
distances than the older calculations. Based on the mod- 
eling literature we pick a value of 120 kpc as an approxi- 
mate distance for the MS-tip (-140°<L M s<-100°). Us- 
ing this constant distance for all the gas in the newly- 
confirmed ~40° portion of the MS-tip gives a mass of 
1.8 x 10 7 M (d/120 kpc) 2 without Wright's Cloud and 
2.8 x 10 7 M (d/120 kpc) 2 with Wright's Cloud. Since 
the mass of the "classical" MS is -5 x 10 8 M (Br05), 
the mass of the new portion of the MS-tip corresponds 
to a ~4% increase in the total mass of the MS. 

As previously mentioned, simulations of the MS indi- 
cate that the distance increases rapidly at the MS-tip 
and, therefore, using a constant distance for the new ex- 
tension is not very realistic. If the distances (as functions 
of Lm s ) from the Kl and GN96 orbits described by |Besla| 
et al. (2007) are used, then the derived mass, without 



WrightVinoud, increases to 5.8 x 10 8 M and 3.0 x 10 7 
M respectively, and 6.3 x 10 8 M and 4.1 x 10 7 M with 
Wright's Cloud. The Kl masses are by far the largest be- 
cause the high- velocity orbit reaches large distances (> 1 
Mpc) for Lms < —125°. Due to ram pressure, the un- 
bound MS gas will slow down, lose angular momentum, 
and fall in its orbit around the MW. Therefore, the ac- 
tual MS distances are likely to be smaller than those from 
the Kl orbit (which is based on the most accurate proper 
motion measurements of the MCs), and a more realistic 
(and conservative) mass estimate is likely to be closer to 
~ 5 x 10 7 M . However, the distances used are poorly 
constrained by observations and the derived masses are 
therefore quite uncertain. 

4. THE ENTIRE MAGELLANIC STREAM 

We combined the MS-tip datacube with the LAB MS 
Gaussians from |Nidever et al.| (2008[) to create images 



of all the detected Magellanic gas. Data were included 
from P03 for small Hi MS cloudets that were not well 
detected in the LAB survey (-95° < L M s < -50°, 
-30° < B M s < -5°). The column density, N Hh of 
the entire ~200° MS+LA system is shown in Figure 
[8^i. The separation of the two main MS filaments at 
Lms^ —80° is clearly visible. The very low column den- 
sity of the MS-tip (L MS < -100°) made it difficult to 
detect in previous surveys, but it is clearly evident in 
Figure [8^i. 

Figure ^jp shows the VLsr _ Lms position- velocity di- 
agram (integrated in Bms f° r —40° < Bms < +40°; 
in units of K deg) for the entire MS. For most of the 
length of the MS (-110°<Lms<-30°) its velocity fol- 
lows a fairly linear gradient (~7.3 km s _1 deg -1 ); this 
changes, however, at the very tip where the velocity lev- 
els out and starts to increase (Lms < — 110°). 

We can now see the WK08 CHVCs (light blue dots in 
Fig. [8^1) and the cloudlets at the far eastern part of the 
MS-tip datacube (near Wright's Cloud) in the context of 
the entire MS. The WK08 CHVCs appear to be an ex- 
tension of the cloudlets near (Lms?L>ms)~( — 45°, +15°), 
just above the main body of the MS, which themselves 
are probably an extension of the Hi filament emanating 
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Fig. 8. — Combined data for all the Magellanic HI gas . (Top) HI column d ensity (log(iV HI ) in units of cm -2 ). The MS-tip datacube 
is used for Lms$ —65° and the LAB MS Gaussians from |Nidever et ah] (2008 ) are used for the rest. The CHVCs from WK08 are shown 
as light-blue dots (not color-coded). (Bottom) Total intensity ot the Magellanic HI i ntegrated along Bms (in units of K deg). The MS 
Gaussians are shown with their true width and not just at their central velocity as in |Nidever et al.f ( |2008| >. The CHVCs from WK08 are 
shown as white dots. The Galactic HI emission is shown as gray contours (at 10, 100, and 10U0 K degj! 
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Fig. 9. — The Magellanic Stream and Clouds in HI (red) with an optical all-sky image (blue, white, and brown; Mellinger 2009) in AitofT 
projection with the direction to the Galactic center at the center. 



diagonally from the center of the Magellanic Intercloud 
Region (ICR) at (Lms 5 ^ms)~( — 9°,— 2°). It appears as 
though all of these cloudlets are connected into one large 
filamentary distribution that is parallel to, but above 
in Bms •> the main body of the MS. This "parallel fila- 
ment" has lower column density and is much clumpier 
than the main MS, and does not appear to continue 
past Wright's Cloud near L M s~ -120°. WK08 pointed 
out that this parallel filament concides with a secondary 
stream predicted b y the numerical simulations by |Gar-| 
diner & Noguchi (1996), although it is less apparent 
or non-existent m more recent models of the MS (e.g., 
Yoshizawa fc Noguchi | [20031 |Connors etliL] [20061 iMalT 
tropietro et al.||2005|). The existence of MS deb: 



2005) 
the MS 



bris on 

Bms < 0°)? which is not seen 



the ot her side o 

in the Gardiner fc Noguchi] ([T996 ) simulations, suggests 
that the MS is generally much broader once low column 
densities are included. Plausibly the process that is cre- 
ating the main body of the MS is also liberating small Hi 
cloudlets that are being spread over a large area of the 
sky. It is therefore possible that even more CHVCs in 
this part of the sky (beyond those discovered by WK08) 
are part of this Magellanic Stream "wake" . 

Figure [9] shows the entire MS with an optical all-sky 
image of the Milky Way and the Magellanic Clouds. The 
large extent of the MS on the sky and its orientation with 
respect to the MW is revealed in this image, as well as 
the fact that the observed MS-tip is close to crossing the 
MW disk plane. 

Figure [lO] shows the column density integrated along 
Bms f° r the entire MS+LA system. There is a strong 
column density gr adient along the MS as has been previ- 
ously noted (Mat hewson et al.]|1974[ P03); this gradient 
is well fitted by JV H i = 53 x lO 21 exp(L MS /19.3°) cm" 2 . 
The LA, on the other hand, has a fairly constant column 
density along its ~60° span. The lower panel shows var- 
ious individual components of the MS and LA: LA I— III, 



and the eastern (S0+S1) and western (S2-S4) portions 
of the MS. 

5. DISCUSSION 



Besla et al. (2007) presented new orbits for the MCs 
that place them much farther from the MW in the past 
than previously thought. If the model orbits are correct, 
then it is difficult to understand how ram pressure or 
tidal forces created a long MS since these forces require 
the MCs to be fairly close to the MW for quite some time 
to be effective. Some other mechanisms, for example, 
dynamical processes associated with star formation, are 
probably required to help the gas escape the MCs when 
at large distances from the MW. The suspected role of 
star formation in the formation of the MS by Besla et al. 
is in agre ement with the in dependent observational find- 
ings of Nidever et al. (2008) that one of the MS filaments 
originates m a region of the LMC with gaseous outflows 
plaus ibly linked to su pergiant shells and star formation. 
Also, |Bekki fc Chiba| (|2009| ) found that they could not 
reproduce the observed MS distribution in their N-body 
models if the new s pace velocities f or the MCs were used. 
On the ot her hand, [Mastropietro ( 2009| ) (and |Mastropi-| 
etro|2010" ) was able to produce a 120 u -long stream using 
ram pressure and the new, higher-velocity MC orbits. 
However, as shown here, the MS is 40° longer than pre- 
viously established — and probably even longer. It is 
unclear from current simulations if ram pressure and/or 
tidal forces alone can account for a stream of this length. 
Thus, our finding of a longer MS offers an additional chal- 
lenge to tidal and ram pressure models. Future surveys 
of the MS-tip may provide even more stringent require- 
ments on MS models and possibly further constrain the 
MC orbits. 

The deviation of the eastern portion of the MS from the 
equator of the MS coordinate system was already seen 
by P03, but the very coherent deviation of the new SO 
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filament for more than —45° has not been seen until now. 
The tidal MS models by Connors et al. ( |2QQ6[ ) reproduce 
this deviation fairly well ( as we ll as multiple filaments), 
while [Mastropietro et aL fl2005|) show a dev iation in the 
opposite direction and |Mastropietr o (2009) show only a 
small deviation (with one filament). It is not entirely 
cl ear why the deviat ion of the SO filament occurs. 
Besla et ah] (|2007[) show that there is a -10° deviation 



of the LMC/SMC orbital paths (to negative Bms) com- 
pared to the location of the MS . Previous MS models 



Gardiner fc Noguchi 1996| |Connors et al. 



2006; 



Mastropietro et al. 2005J) used a small value for the north- 
ern component of the proper motion (/i/v) or the (closely- 
corresponding) X-component of the Galactocentric ve- 
locity (v x ) in order to bring the MC orbital paths into 
alignment with the MS. However, these assumed theo- 
retical va lues are inconsistent with the observed proper 



motions (van der Marel et al. 2002 Kallivayalil et al. 



2006a|b ) . This disagreement between the orbital paths 
and the observed MS cannot be solved by changing the 
MW mass o r by using an aspherical MW potential ( |Besla| 
et al.||2007| ). The newly confirmed -40° extension of the 
MS, and the large deviation of the SO filament (towards 
positive Bms) ? further exacerbates the problem. 

We investigated the Kl and GN96 LMC orbits from 
Besla et al.| ([2007 ) and found that the velocity inflection 
does not correspond to any physically identifiable part 
of these orbit, (unlike, e.g., Vgsr=0 kni s _1 which cor- 
responds to peri- or apo-galacticon) . The inflection is 
just the point in the orbit at which the satellite is ap- 
proaching the observer with the maximum speed, and 
is then dependent on the position of the observer with 
respect to the orbit. However, at a given observing posi- 
tion, the position of the velocity inflection (and the max- 
imum approach velocity) is sensitive to the initial space 
velocity of the computed orbit. This newly found, dis- 
tinct feature in the otherwise remarkably linear velocity 
structure of the MS should provide a new constraint for 
MS simulations. 

We have shown that the velocity of Wright's Cloud is 
in close agreement with the velocity of the MS gas in the 
same region of the sky. This might be an indication that 
Wright's Cloud is part of the MS. Another interpreta- 
tion is that Wright's Cloud is not gas that escaped from 
the MCs, but rather i s part of a "Magellanic Group" 
of galaxies proposed by |D'Onghia fc Lake] ( |2008[ ) to have 
fallen into the MW together. Since the group is proposed 
to be quite extended, some objects will fall in before oth- 
ers, but the whole of the system will enter the MW with 
similar space velocities. Wright's Cloud, therefore, could 
be an Hi cloud (primordial or previously stripped from 
another galaxy) that is falling into the MW behind the 
MCs. Wright's Cloud should then have a similar space 
velocity to the MCs (since they are part of the same 
group) and will experience similar tidal and ram pres- 
sure forces to the MS that should give rise to a radial 
velocity very similar to that of the MS. 

As can be seen in Figure |8j MS debris at low col- 
umn densities is observed at fairly large angular dis- 
tances (-20-30°) from the main body of the MS and 
effectively makes the MS "wider" . This implies that the 
region of the outer gaseous MW halo that has been influ- 
enced by the MS is quite large. The assocation of many 
CHVCs with the MS (WK08 and this paper) makes it 



N(HI) = 5.9x1 21 exp(L MS /19.3°) 




Fig. 10. — Total "cross-sectional" column density, Nui, (in- 
tegrated in Bms) of all the Magellanic HI gas (same as Fig. [8j 
as a function of Lms n units of deg atoms cm -2 ). There is 
a strong column density gradient along the MS which is well fit 
by N m = 5.9 x 10 21 exp(LMs/19.3°) deg atoms cm" 2 . The LA, 
on the other hand, has a fairly constant column density along its 
~60° span. The lower panel shows various individual components 
of the MS and LA: LAI— III, and the eastern (S0+S1) and western 
(S2-S4) portions of the MS. 

likely that yet other CHVCs have MS origins as well. It 
is conceivable that the new Parkes Galactic All- Sky Sur- 



ve y (GA SS; [McClure-Griffiths et al.||2QQ9| 
nd 



al. 12010) and future surveys, such as the 



G 



Kalberla et 
alactic Aus- 

traIIan~5KA Pathfinder survey (GASKAP; J. M. Dickey 
et al. 2010, in preparation), will find that the MS is even 
broader than seen here. 

discovered that at the head of 



Nidever et al. 



^ ( [2008 ) 

the Stream the two filaments of the MS have sinusoidal 
velocity patterns. One possible explanation is that this 
is an imprint of the LMC rotation curve on escaping gas. 
With this hypothesis Nidever et al. were able to measure 
the drift rate of the MS gas away from the LMC to be 
-49 km s _1 and the age of the -100°-long MS (at 50 
kpc) to be -1.7 Gyr. Now that the MS is -40° longer 
the estimate of the age of the MS must also increase. 
Using the same drift rate the new age estimate for the 
MS is at least —2.5 Gyr. The heliocentric distances of 
Magellanic orbits increase significantly at the tail of the 
MS, which in turn increases our MS age estimate via two 
effects. A larger MS distance increases the linear length 
of the MS and probably decreases the drift rate (since 
the ram pressure and tidal forces are weaker at larger 
distances), both of which increase the age. Therefore, 
the MS is possibly even older than —2.5 Gyr. 

The new MS age estimate of —2.5 Gyr closely coin- 
cides with "bursts" of star formatio n in the integrated 



star f ormation hi stories of the LMC ([H arris fc Zaritsky 



2004[ ) and SMC flSmecker-Hane et aLp002| ) at -2-2.5 



Gyr ago; these nearly simultaneous bursts m each sys 
tern are suggestive of an interaction between the MCs at 
that time. The LM C globular cluster age-gap between 



3 and 13 Gyr (e.g., Da Costa 1991 



Rich et aL||2Q01| |PiattTet al.p 



JGeisler et al.||1997 

indicates a dramatic 
onset of star formation in the LMC around —3 Gyr ago, 
altho ugh such an ag e-gap is not seen in th e SMC clus- 
ters J [Da Costap991[ ). Finally, the orbits by |Besla et al.| 
71 indicate possible encounters of the MCs with each 
3 and 6 Gyr ago. Therefore, a plau- 



( 2007 
other around 0.2, 
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sible scenario is that the MCs had a close encounter at 
—2.5-3 Gyr ago that triggered a burst of star formation 
in the LMC and started the self-propagating supergiant 
shell blowout of the MS and LA . 

According to the jBesla et ah] fl2007|) model, the MCs 
were at a distance of ^600 kpc and well outside the MW 
virial radius (~260 kpc) 2.5 Gyr ago. Therefore, an age of 
~2.5 Gyr for the MS implies that the MS started forming 
before the MCs entered the MW. Ram pressure and tidal 
forces would be ineffective at liberating gas at these large 
distances, and this fact further supports the notion that 
star formation feedback and/or a close encounter of the 
MCs with each other are more likely causes for removing 
the MS gas at this stage in its evolution. 

One point that the ram pres sure simulations, t idal sim- 
ulations, and new LMC orbits (Bes la et al.|2 007) agree on 
is that the distance from the MW at the MS-tip (beyond 
^100°) increases quite rapidly (almost exponentially). 
Therefore, the new length of the MS implies that the gas 
at the MS-tip is much farther away than the "classical" 
MS. The MS gas we are observing at Lms ~ 130° is likely 
at a distance on the order of ~200 kpc and possibly even 
larger. This means that the MS spans a larger range of 
MW distances than previously thought, which makes the 
MS a potentially even more useful probe for constrain- 
ing the total MW mass and the MW potential at large 
distances. 

6. SUMMARY 

We conducted an ~200 deg 2 21-cm survey with the 
Green Bank Telescope in the region where the continuity 
of the end of the "classical" MS with the MS-like emis- 



4. The mass of the newly extended ~40° at the MS- 
tip is -2xl0 7 M (d/120 kpc) 2 (including Wright's 
Cloud increases this by ~50%) and increases the 
total mass of the MS by -4%. 



sion reported by Braun & Thilker (2004) was uncertain. 
Our survey, in co mbination with the Arecibo survey by 
|Stanimirovic et aL] (2008), shows that the MS gas is in 
fact both spatially and kinematically continuous across 
this region. We made a combined MS-tip datacube in- 
corporating our GBT datacube, the S08 Arecibo s urvey , 
the BT04 Westerbork data, and the |Brims efall ( |2005 ) 
Parkes data. This MS-tip datacube allows us to draw 
the following conclusions: 

1. The MS is —140° long, approximately 40° longer 
than previously verified. The MS extends to the 
edge of the datacube, implying that the Stream 
may be even longer. Moreo ver, there is one CHVC 
from |de Heij et ah] ( [2002] ) at L M s~ -163° that 
is consistent with an extrapolation of the eastern- 
most filament of the MS in position and velocity. 

2. The tip of the MS is composed of a multitude of 
forks and filaments. S08 previously identified fila- 
ments S1-S4 in the western region, and we identify 
a new filament, SO, on the eastern side. This fila- 
ment splits from SI near (£ms,^ms)~(— 95°, +3°) 
and deviates from the equator of the MS coordinate 
system for more than ~45° until it reaches the end 
of the datacube. 

3. There is an MS velocity inflection near Lms 1 ^ 
— 120° where the MS velocity levels out and starts 
to increase. 



Five C/HVCs from the catalog by |de Heij et al.| ( |2002 ) 
(numbers 237, 304, 307, 391, 402 in their Table 1) match 
the MS in position and velocity and, therefore, might be 
associated with the MS. Wright's Cloud also matches the 
velocity of the MS and is nearby in position, which makes 
an association with the MS possible. The CHVCs that 
WK080 associated with the MS follow the MS velocity 
trend in our Vlsr^ms diagrams and might belong to a 
"parallel filament" next to the MS that stretches from 
the head of the MS to cloudlets near Wright's Cloud. 

The total column density (integrated in Bms) along the 
MS drops markedly and follows an exponential decline 
with the angular distance from the LMC well-fitted by 
N m = 5.9 x 10 21 exp(L M s/19.3°) cm" 2 . 

An increased length of the MS also increases the age 
estimate of the MS. Using the velocity sinusoidal pat- 
tern of the LMC filament of the MS as a chro nometer 
(as was previously done by |Nidever et al. 2008), we es- 
timate that the age of the ~140 u -long Mb is -2.5 Gyr. 
This timescale coincides with bursts of star formation in 
both the LMC and SMC as well as a possible close en- 
counter between the MCs that could have triggered the 
new era of star formation and the formation of the MS 
via supergiant shell blowout. 

These new observational characteristics of the MS 
should offer further constraints on MS simulations and 
provide data that can be used in future work to elucidate 
more accurate dynamical and structural information for 
the Magellanic System and its interaction with the MW. 

We encourage further Hi observations of the MS-tip to 
track the MS to ever greater lengths until its actual end 
is found. 
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APPENDIX 

DETAILS OF THE GBT DATA REDUCTION 

Baseline Removal and The Standing Wave 

During the GBT observing we became aware of a sinusoidal pattern in velocity in the YY polarization that made 
it difficult to use the standard polynomial baseline fitting routines in GBTIDL. The sinusoidal pattern is a known 
standing wave with a period of ~1.5 MHz that arises from a total pathlength of ~200 m and a double reflection (=LA/8 
focus shifts changing the phase by 180°). The standing wave is highly linearly polarized and only appears in the YY 
polarization (Fis her et al.||2QQ3[ ). 
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Fig. Al. — YY polarization data for the GBT brick centered at (7,6)=(103.0°,— 40.0°) showing our GBT baseline reduction steps and the 
sinusoidal pattern. The three-dimensional datacube (position, position, velocity) is shown in a two-dimensional diagram by plotting the 
two spatial dimensions on one axis as the sequential "integration number", (a) The raw YY data; (b) the YY- fitted 5th-order polynomial 
baseline; (c) the YY- fitted sinusoidal baseline model; (d) the polynomial-baseline-subtracted YY data minus the final XX reduced data (to 
remove real emission lines) showing the sinusoidal pattern in the data; (e) the final reduced YY data with both the polynomial baseline 
and the sinusoidal pattern removed - very little residual baseline is apparent. Each diagram has been smoothed with a 20x20 boxcar filter 
to make the fainter features more visible. The units of the color scale are K. The local Galactic zero-velocity gas and intermediate-velocity 
gas is indicated by the red vertical stripe at Channel~350 and the Magellanic Stream is the vertical streak at Channel~580. The MS 
appears clumpy because the integration pattern scans across the MS and the local velocity emission is continuous because it is pervasive. 



As described in §2.2[ we used our own special-purpose IDL routines to perform the baseline fitting and removal. 
The baseline fitting for the XX polarization was fairly straightforward, with a 5th-order polynomial being fitted to 
the 21-cm spectrum (exluding the Galactic emission) with iterative emission line rejection. Removing the sinusoidal 
pattern from the YY polarization was difficult because it was below the noise level of an individual spectrum (~0.12 
K), and therefore multiple steps were required. First, a 5th-order polynomial was iteratively fitted to the spectrum 
excluding Galactic emission and the emission lines previously excluded for the XX polarization. Next, a cosine was 
fitted to the spectrum with the initial polynomial baseline subtracted. The amplitude and wavelength were held fixed 
(A = 0.02 K, A = 305 km s _1 [191 binned channels]) and only the phase was allowed to vary (with initial guess <p = 
192 km s _1 [120 binned channels]). Emission lines were similarly excluded from the fit. Finally, a polynomial -hcosine 
function was fitted to the spectrum using the previously derived parameters as initial guess, and allowing all parameters 
to vary. The cosine parameters were limited to 0.0 < A < 0.04 K, 280 < A < 330 km s" 1 , and 0.0 < <j>< 400 km s" 1 . 
Galactic emission and emission lines previously excluded for the XX polarization and the initial YY cosine fit were 
excluded during this final fitting process. For some of the bricks the phase of the sinusoidal pattern shifted slowly with 
time, but for others it stayed fairly constant. We experimented with many variants of this method and parameter 
constrain ts u ntil we converged on the procedure outlined above, which is adequate for our purposes. 

Figure |A1| shows an example of the YY data for one brick. The three-dimensional datacube (position, position, 
velocity) is shown in a two-dimensional diagram by collapsing the two spatial dimensions on one axis as the sequential 
"integration number". This makes it easier to visualize the overall features in the data. Each diagram has been 
smoothed with a 20x20 boxcar filter to make the fainter features more visible. 

Our baseline removal process is not perfect and there are some wiggles at the mK level. This is probably due to the 
effects of incomplete rejection of channels contaminated by wings of emission lines slightly skewing the baseline fits. 
However, these baseline ripples do not affect our results for the MS. 

The GBT Hi Ghost 



There is a persistent, but weak, emission line permeating our entire GBT datacube (Fig. |A2| left and Fig. A2 right). 
This line does not appear in the S08 Arecibo datacube or in the BT04 Westerbork datacube that both overlap portions 
of our survey and are both more sensitive than our data. Therefore, it is clear that this is not a real emission line, but a 
spurious artifact unique to the GBT data. The line has a Gaussian shape of ^4=0. 004-0.0 13 K, Vlsr =—243.8 km s _1 , 
and a v ~ 7.7 km s _1 , and appears in both polarizations (with equal strength) as well as in the pre-baseline-subtracted, 
frequency-switched data (both positive and negative images). The suspect line maintains a constant velocity and 
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Fig. A2. — (Left) Average spectrum of our entire GBT datacube. The —243.8 km s _1 spurious emission line is clearly visible. A version 
of the average spectrum scaled by 1/1100 and shifted by —240 km s _1 is overplotted in red and is an exact match to the mysterious line. 
This implies that the —243.8 km s _1 line is in fact a "ghost" spectrum in the GBT data. Some baseline ripples are still apparent in our 
data at the mK level, as are a few missed RFI lines. (Right) Intensities of the GBT data, averaged in Galactic longitude, in units of K. 
The left panel shows the data before removal of the —243.8 km s _1 "ghost" emission line, and the right panel shows the data after the 
ghost line is removed by subtraction of a scaled and shifted version of the datacube from itself. The unsmoothed datacube was used for 
the figure and the final Vlsr - & images were boxcar smoothed with a 20x20 filter to enhance fainter features. 

line-width across our entire datacube. The line is not radio frequency interference (RFI) because it is too wide (~200 
unbinned channels) and because it has a constant Vlsr (RFI should have a constant velocity in the observer rest 
frame). The Gaussian shape of the line suggests that it originates from an astrono mical source. However, it is unlikely 
that this is due to the stray radiation from forward spillover of the GBT secondary ( Lockman fc Condon|2Q Q5 ) because 
of the stability of the line (especially in Vlsr and <j v ) over some 200 deg 2 as well as the large negative velocity. 

The only astronomical Hi emission line that is so stable over such a large area is the local, Galactic zero-velocity 
emission. In fact, the line-shape of the —243.8 km s _1 emission is a very close mat ch to that of the zero-velocity emission 



line. The average spectrum of our entire GBT datacube is shown in Figure A2 left with the —243.8 km s emission 



line clearly revealed. We have overplotted a version of the average spectrum scaled by 1/1100 and shifted by —240 
km s _1 that exactly matches the mysterious emission line. This suggests that the —243.8 km s _1 emission line is in fact 
a "ghost" spectrum in the GBT data. We attempted to remove the ghost by subtracting a scaled and shifted version 
of the datacube from itself. The "before" and "after" images of the GBT data (averaged in Galactic longitude) are 
shown in Figure A2 right; the comparison indicates that the removal of the line was successful, leaving no significant 
residual behind. 

We propose that the ghost in our GBT data can be explained as a sideband with constant offset frequency on a local 
oscillator used in the frequency conversion process that was not filtered out properly. This would cause a reproduction 
of the input spectrum at a much lower amplitude and at the offset frequency, in this case ~1.14 MHz, corresponding 
to 240 km s -1 at 21 cm. It is possible that this ghost is present in other GBT Hi data and we suggest that other GBT 
observers watch out for it. 
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